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We argue that an extension of the Minimal Supersymmetric Standard 
Model that gives rise to viable thermal inflation, and so does not suffer from 
a Polonyi/moduli problem, should contain right-handed neutrinos which 
acquire their masses due to the vacuum expectation value of the flaton that 
drives thermal inflation. This strongly disfavours SO(10) Grand Unified 
Theories. The /i-term of the MSSM should also arise due to the vev of 
the flaton. With the extra assumption that m\ — m'jj u < 0, but of course 
m 2 L — ni 2 Hu + |^| 2 > 0, we show that a complicated Affleck-Dine type of 
baryogenesis employing an LH U D-flat direction can naturally generate 
the baryon asymmetry of the Universe. 



1 Introduction 



Thermal inflation [p], |2|, |3|] provides the most compelling solution to the moduli 
(Polonyi) problem 0, [|, ||. However, for a theory of the early Universe to be 
viable it must be capable of producing a baryon asymmetry |7j 

^~3xKr n (1) 

s 

by the time of nucleosynthesis. Thermal inflation probably dilutes any pre- 
existing baryon asymmetry to negligible amounts, and the final reheat temper- 
ature after thermal inflation (Tf ~ few GeV) is probably too low even for elec- 
troweak baryogenesis. Thus if thermal inflation really is the solution of the moduli 
problem, then it is likely also to be responsible for baryogenesis. 

In Section |2| we explain why the flaton that gives rise to thermal inflation 
probably also generates the masses of right-handed neutrinos as well as the \i- 
term of the Minimal Supersymmetric Standard Model (MSSM). We also note 
the various ways in which a potential domain wall problem can be avoided. In 
Section ||] we describe how a somewhat complicated Affleck-Dine type mechanism 
can naturally generate the required baryon asymmetry after thermal inflation. In 
Section £| we give our conclusions. 

2 Thermal inflation, right-handed neutrinos, and 
the //-term 

2.1 Thermal inflation and right-handed neutrinos 

The superpotential of the MSSM is § Q 

W^mssm = XtQH u t + X b QH d b + X T LH d r + fi H H u H d (2) 

Thermal inflation |IJ requires that there is in addition at least one flaton with 
vacuum expectation value |0| = M in the range 10 10 GeV < M < 10 12 GeV, 
the lower bound coming from the requirement that thermal inflation sufficiently 
dilutes the moduli, and the upper bound from the requirement that the final 
reheat temperature after thermal inflation Tf be high enough to thermalise the 
lightest supersymmetric particles (LSP's) produced in the flaton's decay and so 
avoid an excess of LSP's. 

Also, in order for <fi to be held sufficiently strongly^ at <fi = by the finite 
temperature during thermal inflation, <fr must have unsuppressed couplings to at 



1 Here, and throughout most of this paper, all indices (the usual gauge and generation indices, 
as well as any singlet indices) have been suppressed. We will for the most part be focusing on 
the third generation as is suggested by our notation. 

2 Note the strong dependence on Tq in Eq. (35) of Ref. 
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least one other field, say if>, that is light when = 0. We therefore either require a 
term X^cfnf) 2 / 2 with |A^| ~ 1 in the superpotential, or require to spontaneously 
break a continuous gauge symmetry with gauge coupling g^ ~ 1 (with ip being 
the gauge field in this case). One reason to prefer the Yukawa coupling over the 
gauge coupling is that the renormalisation group [§] effect of the Yukawa coupling 
would be to drive the soft supersymmetry breaking mass squared of negative at 
low energies as is required for a flaton, while the gauge coupling would have the 
opposite effect. Another reason is that the gauge symmetry has no independent 
motivation, while, as we shall see, the Yukawa coupling is very well motivated. 
We will therefore focus on the case of the Yukawa coupling. 

After acquires its vacuum expectation value M, if) will acquire a mass 
|A^| M ~ 10 10 GeV to 10 12 GeV and so is not a MSSM field. In order for to be 
coupled to the thermal bath, which is presumably composed of MSSM fields, we 
therefore require if) to couple to the MSSM. We therefore?] require at least oneQ of 
the terms LH u ip or if)H u Hd in the superpotential, as these are the only possible 
renormalisable couplings of a singlet to the MSSM. 

The former is the standard coupling \ v LH u v of a right-handed neutrino if) — v 
to the MSSM, and furthermore v automatically acquires a mass M v = |A^| M ~ 
10 10 GeV to 10 12 GeV in the right range for the seesaw mechanism |TIJ to generate 
a left-handed tau neutrino mass 

...j, r „ (174 GeV) 2 sin 2 /? 

m "L - — - i\ Hz y 6 ) 

—) (4) 
0.5 J K J 

suitable for the mixed dark matter scenario for the formation of the large scale 
structure of the Universe |TT| . 




The latter coupling if)H u Hd does not have any good independent motivation, 
apart from the fact that it is possible. We will therefore focus on the former case. 
Later though, we shall make use of this possibility for coupling to the MSSM in 
a different context. 



3 Assuming ip is not charged with respect to the MSSM continuous gauge symmetries as this 
would, in general, destroy susy GUT gauge coupling unification. However, if ip is a complete 
SU(5) multiplet the unification of the gauge couplings will be unaffected. For example, one 
could have A^i/j 2 = A^055. Alternatively, appropriate choices of representations could shift 
the unification scale to the string scale ||. 

4 If we have both then we need two different "0's, one charged and one neutral under R-parity. 
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2.2 The final reheat temperature after thermal inflation 
and the /i-term 

Now that we have a more precise picture of the couplings of the flaton, 



so far 



Wi 



(5) 



we can hope to make a more precise estimate of the decay rate of the flaton and 
so the final reheat temperature after thermal inflation Tf. 
We first note that the effective superpotential coupling 



W m 



(X U LH U ) 



(6) 



obtained by integrating out u, i.e. eliminating v via the constraint dWjdv = 0, 
will give a decay rate of order T ~ m 5 /M 4 which is negligible. 

The dependence of the low energy renormalised coupling constants will give 
larger decay rates. To estimate these we first need to know the contributions of 



the right-handed neutrinos to the renormalisation group equations. They are [12 



at 

16vr 2 4-A T 
dt 

i« 2 d 

167T —Jl H 

dt 

16n 2 — rn? 
dt L 

i a 2 d 2 

Io7T —rrij 



|AJ^ + ... 



I At/ 1 A T + . 



|Aj/| + . . . 



2 d 

1671 Jt » HH " Hd 



2|A, 
2|A t 
2|A t 
2|A t 
2IA, 



(m| + m 2 + + \ A LHuV \ 2 ) + 



Alh u v + ■ ■ ■ 



AlH u v + • • • 

^4l_h" U !/ + • • • 



(7) 
(8) 
(9) 
(10) 

(11) 
(12) 

(13) 

(14) 



where t is the logarithm of the renormalistion scale, the m's are the soft supersym- 
metry breaking masses, the A's are the soft supersymmetry breaking parameters 
in the terms in the scalar potential of the form AW + c.c, the m's and the 
magnitudes of the A's are of the order of the soft supersymmetry breaking scale 
m s ~ 10 2 to 10 3 GeV, and the . . . stand for other terms independent of the right- 
handed neutrinos. \<p\ sets the threshold for the right-handed neutrinos, and so 
writing \(f>\ = M + 5<f> r /y2 we get the effective couplings 

|2 



w. 



A„ 



cff 



16^2 tt 2 M 



(X t QH u t + X T LH d r + n H H u H d ) 



(15) 
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and 



|A 12 



^softeff = 8v /2^2 M + m l + m H u + \ A LH u v\ 2 ) (}L\ 2 + \H U \ 



IAJ 2 



+ o at ,^ A LHuu (XtQH u t + X T LH d r + fi H H u H d ) 5<j) r (16) 
ov 2 tc z M 



From these couplings we estimate the total decay rate to be 

I \ i4 3 

10 4 M 2 

This would give a final reheat temperature after thermal inflation of 



(17) 



g* J T*M^ (1 



'3 x 10 12 GeV\ / m 



~ 10 MeV ^— ^— (19) 

V M/\X U \ 2 J V 300 GeV 7 1 ; 

where the first bracket is constrained to be of order one, or perhaps less, by 
Eq. (f|). However, in order not to over-produce LSP's we require |L3) 

m \ 1.5 to 2 



T f > 1 GeV -= — (20) 

V 300 GeV/ v ; 

Thus our model seems to be in trouble unless we can add some extra coupling 
that gives a stronger decay rate. The only possibility is to couple (f> to H u H d in 
the superpotential. A term (f)H u H d would require a very small coupling constant 
to avoid generating too large a /i-term, but a termQ 

^decay = (21) 

is not only allowed but could naturally generate a /i-term 

^ - -j^r (22) 

of the required size |14" |. For example, for M = 10 11 GeV and |A M | = 0.1 we get 
l/i^l = 400 GeV. From now on we will assume that the //-term is generated in this 
way so that fi H = 0, or at least \fi H \ < |/^|, and |/^| ~ m s . 

Writing <p = M + 5<p, where |.M| = M and /i^ = A M X 2 /M P1 , we get the 
relatively unsuppressed couplings 

£decay = 2^<pH u H d j^ — 2 |/i^| 2 (j-ff«| 2 + |-^d| 2 ) 

- 2^ (^Qt#d + + VM U + + c.c. (23) 



5 A term X^HuHd/M^^ would be an alternative if M > 3 x 10 11 GeV. We assume the 
displayed case for simplicity. 
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where a tilde denotes the fermionic component of the superfield, a bar denotes 
the hermitian conjugate, and is the soft super symmetry breaking parameter 
in V so ft = A^n^H u H d + c.c. and so \A^\ ~ m s . The decay rate is then estimated 



to be 

„ mz, 



10 2 M 2 ^ 
where we have roughly assumed ~ ~ \A^\ ~ m B . We therefore get a 
reheat temperature 

'(™)(30^V) § 

which is sufficiently high. 

We thus expect the flaton to have the following couplings to the MSSM Q 

^couplings = X v LH u v + ^A^z/ 2 + H^Hd ^ 

z ivipi 

2.3 The flaton potential and domain walls 

In this section we consider the self-couplings of the flaton. The flaton (or, in the 
case of a multi-component flaton, at least one component of the flaton) should 
have a negative soft supersymmetry breaking mass squared — m 2 |0| 2 to drive it 
away from = after thermal inflation. It will also need a term to stabilise 
its potential at its vacuum expectation value \<p\ — M ~ 10 10 to 10 12 GeV. The 
simplest possibility is[] 

and this is what we shall assume. In the case of a multi-component flaton this 
could be interpreted as for example AAf0i/4M P1 , \ M cj)i<^/ My>\ or a sum of such 
terms. See Ref. for an explicit multi-component example. Here, for simplicity, 



we will focus on the case of a single component flaton, though it should be born in 
mind that a multi-component flaton might be preferable from the model building 
point of view. 

We then get the following scalar potential 




V(<P) = V - ml H 2 + ( ——77—^— + c.c. ) + (28) 



pi 



6 As mentioned before, A^ 3 H u Hd/ 'Mq VT is a possible alternative to X^(f> 2 H u Hd/Mp\. 

7 W vcv = \m4> 5 /Mgut wou ld be an alternative if M > 3 x 10 11 GeV. Again we assume the 
displayed case for simplicity. One might even be able to use the renormalisation group running 
of rrifj, to stabilise the potential but then one would not automatically get a value for M in the 
correct range. 
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where m^ ~ \Am\ ~ m s . This potential has four degenerate minima with \<f)\ = M, 
where 



m2 = 2m M P1 / | A M | 



3 A 



M\ 



\A M \ 2 | 3 
4 



(29) 



For example, for = \Am\ = 300 GeV and |Am| = 0.1, we get M = 10 11 GeV. 
The eigenvalues of the mass squared matrix at the minima are 



lQmJA 



m. 



4>\-n-M\ 



\A 



M I 



+ 



m<t> \ 



\A M \ 2 | 3 
ml A 



(30) 



and 



m A 



Ami + mjL 



Ami 



1 + 



A\A 



M 



3m rf 



1-4 



M 



+ 



m^ \ 



\A M \ 2 | 3 
ml A 



(31) 



and the flatino mass squared is 



\A 



M 



+ 



m 



Requiring zero cosmological constant at the minima gives 



\A M ? 3 
2 + 4 



(32) 



V = -m\M 2 



1 + 



2\A 



Ml 



\A 



3m d 



m^ 



+ 



\ 




(33) 



With four degenerate minima we clearly have to worry about a potential do- 
main wall problem. The simplest way to eliminate the domain walls is to add a 
small term which breaks the degeneracy of the vacua, the difference in pressure 
exerted on the walls causing the domains with greater vacuum energy to collapse. 
They collapse before the domain walls come to dominate the energy density if 
the difference in vacuum energies e satisfies e > <r 2 /Mp 1 where a is the energy per 
unit area of the domain walls [IB|. For flaton domain walls a 
require 



m s M and so we 



e > 



m 2 s M 4 



(34) 



Therefore a term in the superpotential of the form 



m s M 4 - n (p n 

Mi 



Wwalls > 



(35) 



with n odd would be sufficient to eliminate the domain walls. A term with 
n = 2 mod 4, for example W^aiis ~ 6 /^pi> would reduce the Z 4 domain walls 
to Z 2 domain walls. Note that W^is can be extremely small, and hence have 
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a negligible effect on the dynamics to be discussed in the next section, but still 
solve the domain wall problem. 

Another way to avoid a domain wall problem is to gauge the discrete symmetry 
so that there is really only one vacuum. However, non-trivial anomaly cancellation 
conditions must be satisfied [|I7|. In the case of a single-component flaton with 
the superpotential of Eq. fl36l) , and no extra light SU(3) multiplets, the mixed 
discrete-SU(3) anomaly cancellation condition requires <p 2 to be neutral under 
any unbrokenP] anomaly free discrete gauge symmetry. We therefore cannot use a 
discrete gauge symmetry to remove the Z4 domain walls. However, the anomaly 
free Z 4 subgroup of XJ(1)b-l+4Y, under which has charge 2, can be used to 
gauge away the Z 2 domain walls left by W wa u s ~ 6 /M|j above. Furthermore, 
this symmetry is broken down to the standard matter parity of the MSSM (which 
is equivalent to the R-parity of the MSSM) by the vacuum expectation value of 

In the case of a multi- component flaton the discrete symmetry can be extended 
to a continuous symmetry which may or may not be gauged. In the case of a 



continuous global symmetry, for example the Peccei-Quinn symmetry of Ref. [|T5 
the Goldstone bosons may prove troublesome JTg| . One also has more freedom 
to satisfy the anomaly cancellation conditions in the case of a multi-component 
flaton. 



2.4 Summary 

An extension of the MSSM that gives rise to viable thermal inflation, and so 
does not suffer from a moduli problem, should have the following terms in its 
superpotential 



W ti = X t QHj+X b QH d b+X T LH d T+X u LH^+l\^u 2 + X ^ uH U ^^ (36) 



3 LH U Affleck-Dine baryogenesis after thermal 
inflation 

To orientate the reader we will first sketch the basic idea we have in mind before 
plunging into the details. 

The .D-flat direction parametrised by LH U provides an ideal Affleck-Dine field 



1^, [20|, |5[. In order for it to behave as an Affleck-Dine field we must first get 
it away from zero. We therefore require its mass squared at the end of thermal 
inflation, [m 2 L — m 2 Hu + \^h\) /2, to be negative. It is simplest to assume that 



it rolls away from zero before cf> does. However, when = the right-handed 



R-symmetries are broken down to Z2 by hidden sector supersymmetry breaking. 
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neutrinos are light, and so LH U is not F-flat (it has a quartic term in its potential 
coming from the superpotential) . LH U will thus be stabilised at a modest value. 

Next will roll away from zero. The right-handed neutrinos become heavy 
and so can be integrated out leaving the effective seesaw coupling W seesaM given 
in Eq. (||). This is now the term that stabilises the LH U direction and we see that 
it gets smaller, and so the LH U direction gets flatter, as gets larger. Thus, as 
rolls away from zero, LH U will roll further away from zero. Furthermore, the soft 
supersymmetry breaking term derived from W^eesaw will correlate the phases of 
and LH U . 

When becomes sufficiently large, it will start to feel the basin of attraction 
of one of the minima of its potential, and so will start curving in towards that 
minimum, i.e. its phase will be roughly determined modulo tt/2. The phase of 
LH U will then be roughly determined modulo 7r/4 by the soft supersymmetry 
breaking term derived from W seesaw . 

When |0| becomes of order M, a cross term from the supersymmetric part of 
the potential becomes significant and changes the correlation between the phases 
of and LH U , and so gives the phase of LH U a kick. The direction of the kick is 
determined by the parameters in the lagrangian (this is our CP violation) and so 
gives a non-zero net contribution when averaged over different spatial locations, 
unlike the rest of the angular momentum that is flying around. Furthermore, as 
this is happening Wdecay (see Eq. fl2T|)) starts to give a significant contribution to 
the mass of H u and hence LH U . For m 2 L — m 2 Hu + \hh + H<p\ 2 > this gives the 
LH U direction an overall positive mass squared (as it must because LH U has a 
positive mass squared in the true vacuum), and so sends LH U spiralling back in 
towards zero. 

The effective friction on the motion of and LH U coming from the Hubble 
expansion is negligible. However, the effective mass squareds of both and LH U 
have been changing sign during the above dynamics and so one would expect 
them both to decay via broad parametric resonance |21 ]. This will lead to ap- 
proximately critical damping, and so it seems reasonable to expect that both 
and LH U will be trapped near their vacuum expectation values essentially imme- 
diately after the dynamics described above has occurred. Once they are trapped, 
parametric resonance becomes less efficient because the mass squareds are now 
always positive. LHJs potential near LH U = conserves angular momentum, 
or in other words lepton number, and so LHJs newly acquired lepton number is 
conserved. 

The dynamics outlined above is illustrated in Figure 1. 

The decay of the LH U Affleck-Dine condensate will generate enough partial 
reheating to restore the electroweak symmetry, and so its lepton number can be 
converted to baryon number by the usual electroweak effects [p2| |. Note that the 
energy density is still dominated by the flaton and the reheating in the Affleck- 
Dine sector has a negligible effect on the now decoupled flaton. Finally, after the 
temperature has dropped to a few GeV, the flaton decay will complete, releasing 
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substantial entropy. 



3.1 Estimating the baryon asymmetry 

Our basic model is 



Wti = XtQHj+\ b QH d b+\ T LH d T+KLH u u+^ ( l ) u 2 + X ^ uHd + ^^ (37) 



Mi 



pi 



4Mi 



pi 



The squark fields have no linear terms in their potential and have positive mass 
squareds. They will therefore be held at zero apart from thermal fluctuations, 
and so can be ignored apart from their contribution to the finite temperature 
effective potential. The zero temperature potential for the other fields is 



V 



\X T LH d \ + \\ V LH U + A^0z/| + \\ T H d T + \ u H u v\ + 



\^ 2 H d 

\ v Lv + 



M 



pi 



+ 



\ t Lt + 



M 



pi 



+ 



1 x 2 . n X^H u H d A M </> 3 

-A^Z/ + Z — r 



M 



PI 



M 



pi 



+ .D-terms 



( a \ t tt , a \ t u I a \ / 2 i A^ii4> 2 H u H d A m Xm<P 
+ A T \ T LH d T + A v \ v LE u v + A^A^^z/ + - H 



M 



i 2 I U , 2 I l-a i 2 I r I 



PI 



Mpi 

2 I jl|2 



C.C. 



m 2 M 



(38) 



where the m's and the magnitudes of the A's are of order m s . We assume 



LH U 



(39) 



so that the D-flat direction parametrised by LH U is also unstable, in addition 
to the flaton 0. Note that after <fi acquires its vacuum expectation value M ~ 
m s Mpi/\\M |, it will give an extra contribution |A M | 2 M 4 /Mpj to HJs mass squared. 



This will be of order m? if |A,, 



Am I- We assume 



m 



LH U 



=M 



1 

2 



m, 



m 



tin 



+ 



|A M | M 4 
M^ 



> 



(40) 



so that the direction is stable in the true vacuum. 

A rigorous study of the dynamics of this model is beyond the scope of this 
paper. Instead we will make some simplifying assumptions in order to illustrate 
how the Affleck-Dine mechanism might be implemented after thermal inflation 
and to crudely estimate the resultant baryon asymmetry. 

We assume that all fields are initially held at zero by the finite temperature 
during thermal inflation. We assume that the LH U direction rolls away from 
zero first. It will be quickly stabilised by the term \\ U LH U \ 2 at a value \LH U \ ~ 
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m 2 LH ^/\X v \ 2 . The term A v X v LH u v + c.c. then causes v to roll away from zero. 
Then the term X u LH u X^u(j) + c.c. causes to roll away from zero[]in the direction 



argi 



arg 



X$A V (X U LH V 



(41) 



For simplicity we will assume that r and remain at zero, or at least that any 
expectation values they acquire can be neglected. With this assumption, once 
and LH U escape beyond the reach of the temperature, their dynamics will be 
governed by the zero temperature potential 



V = \X U LH U + A«i0z/| 2 + \X u H u i>\ 2 + \X u Lu\ 2 + 



X^ 2 H U 



Mi 



pi 



+ D-terms + A^X u LH u u + A^X^v 2 + 



AmXm4> 4 



M 



+ 
+ c.c. 



1, 2 . A A /0 3 
-AaV H 

2 Mpi 



pi 



i 2 I |2 | 2 I r |2 2 tt i 



(42) 



As |0| increases, v will quickly acquire a large mass ~ |A^0|, and so will be 
constrained to the minimum of its potential 



v 



X U LH U 



(43) 



The effective potential then becomes 



V = (\X V H U \ 2 + \X U L\ 2 ) 



X U LH U 



X^ 2 H U 



M, 



pi 



+ 



A M 3 



M 



pi 



+ D-terms 



+ 



As - A v + 



Am0 3 \ (X u LH u ) A m X M (, 



+ 



M 



+ c.c. 



pi 



i 2 I t |2 2 I tt |2 2 



(44) 



We assume the D-terms constrain L and H u to the D-flat direction LH U . Then 
writing LH U = ip 2 /2 we get 



V 



A^ 3 



2A^0 
+ 



+ 



A M 2 V 



Mi 



pi 



+ 



Am0 3 



Mpi 

,2„/,4 



- m„ 



Ad, - A„ + 



A M 3 \ A 2 V , A M A M c 



20 



+ 



M 



I 2 - ml 



+ c.c. 



pi 



(45) 



9 One might imagine that our Affleck-Dine type mechanism could also be implemented using 
say the right-handed electron sneutrino, which could plausibly have a small quartic coupling 
A^e, instead of LH U . However, unlike LH U , if it was unstable it would roll away from zero at 
some early time because all its couplings would be small. The term A^X^v 2 4> + c.c. would then 
cause 4> to roll away from zero causing a premature end to thermal inflation. 
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\ — m 2 L )/2. To make this potential more transparent, we make 



where m\ = [m 
the following change of variables 

V = m 2 M 2 V , = Mcp , ip = Mip , = m s a0 , = m s a^ (46) 

m s Mpi 



, , , x • - . . , x m s M Pl A M 

A v - A 4> = m s a u , A M = m B a M , A M = p- , \ M = , , 9 Pm , 77 = — 7 

A£ m s 

(47) 

where the a's and the magnitudes of the a's and /3's are of order one and we 
assume |7| <C 1. We then get 



V 



~ 2 



+ |/?Af | 



+ 



(-4 + 


</> 2 




270 



+ 5 Ifl., 







• 




2 










+ 



( 



oimPm <P ~ 



_ ~4" 
«m/?A/0 



V 



«M0 



470 



- + c.c. 



(4£ 



When <C 1, "0's potential is stabilised at 

~ 2 

1p ~ |7| 

while its phase is coupled to that of by the termQ 



(49) 



ocm- 



_ ~4 
«M/?M0 



\ 



a M (p 



470 



- + c.c. 



(50) 



the second term in the brackets being negligible at this stage. 

When 3 > 7, the potential for the phase of will be dominated by the term 



«M/?M0 4 + C.C. 



(51) 



and so in some sense we can regard the phase of as being determined modulo 
7r/2. Put in a different way, will be pulled towards one of the minima of its 
potential and so its phase will be strongly biased towards 



«m/5m0 4 = - o> m (3m4> 4 



(52) 



The phase of ip is then determined modulo ir/8 by the term in Eq. (|50f ). 

When becomes of order one, two things happen. First, the second term in 
the brackets in Eq. (|50|) becomes of order one and gives the phase of ip a kick in 



10 The correlation induced by this term is different from that of Eq. ( p~[| ) and so the phase of 
ijj will get a kick in the direction sin(arg A$ — arg-Aj,) while the phase of <f> will get a kick in the 
opposite direction. This may contribute to the net lepton number generated, in addition to the 
similar effect to be described below. 
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the direction sin(arga I y — arg«M)- Note that even before this ip will have had 
some angular momentum about ■0 = 0, but it averages out to zero in the Universe 
as a whole, as is shown in Figure 2(a). This new contribution has a direction 
determined by the parameters of the lagrangian and so will give a non-zero net 
contribution, as is shown in Figure 2(b). Put another way, the difference in phase 
between a v and au is our source of CP violation. Secondly, the last term in the 
brackets in 



-tty + 



270 







1 




2 










(53) 



becomes of order one giving -0 a net positive mass squared and so causing it to 
spiral back in to -0 = 0. 

Assuming the expected broad parametric resonance |2T] provides enough damp- 
ing, both and ip should then become trapped near their vacuum expectation 
values, after which the parametric resonance becomes less efficient. ?/>'s potential 
near -0 = conserves angular momentum, or in other words lepton number, and 
so ?/>'s lepton number is conserved. 

The Affleck-Dine condensate ip will decay well before the Hubble expansion 
reduces its amplitude to the electroweak scale, and so will release enough thermal 
energy to restore the electroweak symmetry. The lepton asymmetry will then be 
converted into a baryon asymmetry 



n B 
s 



1 (ul 
3 V s 



(54) 



by the usual electroweak effects Finally, after the temperature has 

dropped to a few GeV, the flaton decay will complete, releasing substantial en- 
tropy. 

The baryon asymmetry generated in this way is roughly estimated to be 
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T { n L 6T f n^ 



9\X 
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0(100 GeV) 2 T f 

77V M 2 




(55) 
(56) 



where we have used Eqs. (f49|) and (|3|). 6* is defined by this equation and will depend 
on the phase difference between a v and olm as well as the detailed dynamics. 
As discussed in Section 0, we expect Tf ~ 1 to 10 GeV and M ~ 10 10 to 



10 GeV. Neutrino phenomenology |LT| suggests m Vr 



to 10 3 eV and m VeL ^^^ Mi - Therefore, in order for Eq. 
asymmetry of Eq. ([!]), we require 9 to be roughly^] 

6 T ~ 10" 4 to 10 



5 eV, m VfML ~ 



io- 



6j) to give the baryon 



(57) 



> 1 corresponds to the scenario being unviable. 
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# e < # M ~ 1(T 7 to 1(T 2 (58) 

depending on which generations make up the Affleck-Dine LH U direction. The 
eventual measurement of the Higgs and slepton masses should help to determine 
which of these ranges is the appropriate one (or rule out the whole scenario), and 
a measurement of m„ eI would narrow the uncertainty in 9 e . 

4 Conclusions 

Right-handed neutrinos should acquire their masses due to the vacuum expecta- 
tion value of the flaton that gives rise to thermal inflation, not some composite 
GUT operator. This will have important implications for GUT model building. In 
particular, SO(10) GUT's are strongly disfavoured because the flaton would have 
to be in a 126 representation which is difficult to derive from superstrings and one 
would have a flaton-125 splitting problem in addition to the usual doublet-triplet 
splitting problem. 

The //-term of the MSSM should also be generated by the vev of the flaton. 
Our Affleck-Dine type mechanism generates a baryon asymmetry which is 
roughly estimated to be 

where 9 is the lepton asymmetry per Affleck-Dine particle. 9 depends on the dif- 
ference in phase between the soft supersymmetry breaking parameters of PU see saw 
and W vev (c.f. Eqs. (|6|) and fl2"T|)), as well as the detailed dynamics. 
We also make the following prediction 

m i < m # (60) 

modulo renormalisation effects, where — m 2 H is the soft supersymmetry breaking 
mass squared of H u , and m 2 L is the soft supersymmetry breaking mass squared of 
a lepton doublet. 
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Figure captions 



Figure 1 

Numerical simulation to illustrate the dynamics of <f> and LH U = if) 2 /2 after 
thermal inflation. The potential of Eq. ( fi8| ) was used with the parameters a<p = 
= a u = (3m = 1) \o£m\ — — 2, 7 = 10~ 3 and arga M = 3.1. The initial 
conditions were \<f>\ = = 10~ 3 , arg</> = 1 and aigip = 0.25. A friction term Tip 
with r = 0.75 was added to the equation of motion of to crudely simulate the 
effects of parametric resonance. A friction term was not added for ip because it 
would obscure the total lepton number generated which in reality is contained in 
both the homogeneous ip field and its decay products (such as the inhomogeneous 
ip modes produced by parametric resonance). 



Figure 2 

Numerical simulation to show the non-zero net lepton number generated. The 
same parameters as in Figure 1 were used except for the following. Motivated by 
Eq. (|4T|), the initial phase of ip was taken to be random while the initial phase of <p 
was taken to be given by arg(0) = 4 arg(^) + C . The lepton number produced, as 
measured by &rgip(t = 100) — arg = 0), is plotted against the initial phase of 
if). The dotted line gives the average value. The plots correspond to the following 
values of the parameters (a) arga^ = ^ (which is CP conserving) and C — 0, 
(b) arg«M = 3.1 and C — 0. 
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